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form of the claims for examination. It is respectfully requested that examination be 
performed on claims 1 1-38. 
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INSTITUT FRANQAIS DU PETROLE 

METHOD FOR DETERMINING A SEISMIC WAVE VELOCITY MODEL 
IN A HETEROGENEOUS UNDERGROUND FORMATION 

Inventor : Karine BROTO 

10 ABSTRACT 

FMl Mothod A method f or determining a velocity model of seismic waves 
picked up by seismic receivers coupled with an underground formation, in response 
to the emission of seismic waves in the subsoil by a seismic source, after reflection 
on geologic interfaces of said -the f ormation, from multi-offset records of these waves 

15 is disclosed . f[- ]]From picking ©ivthe seismic records of one or more events for 
different source-pickup pairs associated with the weakest offset of a seismic 
acquisition device, the geometry of the events con si d e r e d is predicted for the total 
range of offsets of the acquisition device and for different velocity values taken 
around a reference layer velocity, by means of kinematic inversion and modelling 

20 techniques. Then, for different lateral positions of the acquisition device, the 
predicted travel Jime curve best approaching the seismic event considered is 
selected. The traveljimes thus obtained are then supplied to a prestack inversion 
method such as a traveUime tomography. 

App li cat i ons : s e i sm i c s ub s oi l i mag e ry for e xamp le . 



BACKGROUND OF THE INVENTION 

F I ELD OF THE I NVENT I ON Field of the Invention 

The present Invention relates to a method allowing to d e torm i ne determination 
of a velocity model of seismic waves picked up by receivers coupled with an 
5 underground formation, from multi-offset records of these waves. 

Th e m e thods a ll ows eas ie r acc ess to th e pr es t a ck k i n e matic i nformat i on 
assoc i ated w i th the ovonts conta i n e d i n th e s e ism i c r e cords, wh i ch i s comb i n e d w i th 
a prostack kin e matic i nv e rsion t e chniqu e for i mag e ry of th e g e o l og i c i nt e rfac e s of tho 
subso i l, 

10 What is r e f e rr e d to as "pr e stack k i n e matic informat i on" i s th e trav el t i mes 

as s ociat e d w i th th e r e fl e ct i ons r e cord e d by sourc e p i ckup pa i rs locat e d at var i ab le 
distanc e s from on e anoth e r W e contrast " e asi e r acc e ss" w i th manual p i ck i ng of the 
soicm i c ovonts i n tho mu l ti - offs e t co ll ect i ons. A "pr e stack kin e mat i c i nversion 
m e thod" i s und e rstood to b e a m e thod a l low i ng, from th e k i n e mat i c i nformation 

IS e xtract e d from the se i sm i c r e cords (and not from an approximation of this kin e matic 
i nformat i on), to f i nd the geometry of the assoc i ated r e fl e ctor s and th e l ay e r 
ve l ocitios. 

BACKGROUND OF THE I NVENTION 

DESCRIPTION OF THE PRIOR ART 

20 Seismic reflection surveys are widely used in petroleum exploration, notably to 

produce images of the subsurface from the information contained in the waves 
propagated and reflected on -from t he geologic discontinuities of the subsoil. 
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More precisely, imagery methods use the kinematic information associated 
with the major seismic reflections (irer -that is t he travel.times of the waves reflected 
on the main discontinuities of the subsoil) to determine a velocity macromodel of the 
subsoil, which will be used to convert the temporal seismic records to a depth image 
5 of the subsurface. 

Access to the kinematic information necessary for determination of the 
velocity model requires interpretation of the seismic events in the seismic multi-offset 
records. Now, seismic multi-offset data is generally characterized by a bad signal-to- 
noise ratio, hence the failure of conventional automatic picking and the necessity of 
10 carrying out a long and costly manual picking of the seismic data. Besides, for 3D 
seismic surveys (currently predominant in relation to 2D surveys), the quantity of 
data to be interpreted is a 4D volume that can reach sizes of the order of one 
hundred gigabytes, or even of the order of a terabyte, which makes the interpretative 
task even longer and tedious. 

15 In order to avoid this stage of seismic multi-offset collection interpretation, 

geophysicists have developed methods based on approximations of the geometry of 
the seismic events in the multi-offset collections. To establish these approximations, 
these methods put fonA^ard hypotheses on the subsoil complexity, hypotheses which 
can relate to the geometry of the geologic discontinuities of the subsoil and to the 

20 layer velocity variations. The method described by Taner and Kohler (1969) can for 
example be mentioned, which assumes a stratified medium con sis ting of horizontal 
plane layers with homogeneous and isotropic layer velocities, as well as the method 
described by Levin (1971) which extends the previous method to sloping layers. 
Other variants have been proposed, but generally the existing methods are based on 



a hyperbolic hypothesis of the traveljime curve in the multi-offset collections. This 
hypothesis is however violated from the moment that the geologic discontinuities of 
the subsoil are no longer plane and/or that the layer velocities exhibit lateral 
variations. 

5 SUMMARY OF THE INVENTION 

W e propo se h e r e The invention is a method for best approaching the travel 
times of the multi-offset collections in cases where the geologic discontinuities of the 
subsoil have any geometry and where the layer velocities are moderately variable 
laterally, while limiting the human time required for interpretation of the seismic 
10 records. The traveljimes obtained are then treated by means of a prestack 

kinematic inversion method a ll ow i ng to tak e taking account of complex kinematics, 
such as travel.time tomography, described for example in the following document: 

[[- ]]Bishop, T. et al., 1985, "Tomographic Determination of Velocity and Depth in 
Laterally Varying Media", Geophysics, 50 N0..6, 903-923. 

15 The method can also apply in cases where the layer velocities are laterally 

variable, but iteratively, the velocity model after updating at a given iteration by 
applying a prestack kinematic inversion method such as the aforementioned 
reflection tomographic method serving as the initial model for the next iteration. 

The method of the invention allows easier access to the prestack kinematic 
20 information associated with the events contained in the seismic records, which is 
combined with a prestack kinematic inversion technioue for imagery of the oeolocic 
interfaces of the subsoil. 
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What is referred to as "prestack kinematic information" is the travel times 
associated with the reflections recorded by source-pickup pairs located at variable 
distances from one another. Contrasting of "easier access" with manual picking of 
the seismic events in the multi-offset collections is performed. A "prestack kinematic 
5 inversion method" is understood to be a method allowing, from the kinematic 

information extracted from the seismic records (and not from an approximation of 
this kinematic information), to find the geometry of the associated reflectors and the 
layer velocities. 



10 SUMMARY OF THE I NVENT I ON 

The method according to the invention allows to d e t e rm i n e determination a 
velocity model of seismic waves picked up by seismic receivers coupled with an 
underground formation, in response to the emission of seismic waves in the subsoil 
by a seismic source, after reflection en-bv geologic interfaces of said -the f ormation. 
15 from multi-offset records of these waves. 

For each seismic event located on-in the records and for each layer delimited 
by said-interfaces, the method comprises at least the following st a ge s steps : 

a) from the prestack seismic records, constructing an iso-offset collection 
from which kinematic information or traveljimes associated with the event are 

20 extractedri 

b) selecting a velocity range around a reference velocity in said -the layer, that 
is sampled with a predetermined intervalTi 
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c) for each velocity sample, applying an inversion technique at affixed 
velocity so as to determine, from the traveljimes extracted from the iso-offset 
collection, the-a_geometry of salct -the interface for the velocity sample concerned in 
order to obtain a series of interface/velocity pairs for said -the e ventv: 

S d) calculating the-kinematic information associated with each 

interface/velocity pair obtained, for source-receiver pairs corresponding to multi- 
offset collections existing in the seismic recordSri 

e) for each interface/velocity pair and for each multi-offset collection selected, 
evaluating the coherence between the multi-offset travel Jimes thus calculated and 

10 the seismic records, and selecting for each multi-offset collection the travel Jime 
curve showing maximum coherence with the seismic recordsri 

f) applying a prestack kinematic inversion method using the multi-offset travel 
times obtained for all the multi-offset collections selected, in order to determine the 
geometry and the velocity of the layer consideredTi and 

IS g) iterating n times (n ^ 0) stag e s steps a) to f) by considering on each 

iteration the velocity model obtained during the previous iteration as the reference 
model to define the reference velocity of the new iteration. 

According to an implementation mode suitable in cases where the velocity 
range selected is not precise enough at the end of either the previous iteration or of 
20 s tages step a) to f), stage -step g ) is carried out n times with n ^ 1 . 

According to another implementation mode suitable in cases where the 
ve l oc i t i es v elocitv d istribution varies greatly laterally and/or in cases where no 
sufficiently precise a priori knowledge of the velocity distribution in the layer 



considered is available, stage -step g) is carried out on offset ranges that are 
increasingly wtde HCireater a s iterations progress. 

According to another implementation mode suitable in cases where the 
ve l oc i t ie s velocity d istribution varies greatly laterally and/or in cases where no 
S sufficiently precise a priori knowledge of the velocity distribution in the layer 

considered is available, staoe -step g) is carried out on multi-offset collection grids 
that are increasingly finer as iterations progress. 

According to an implementation mode suitable in cases where the s ought 
interface geometry which is sought g enerates triplications, stage -step e ) is carried 
10 out by considering the triplication branches in the multi-offset collections 

independently of one another. For example. R ay -rav t racing and inversion tools 
allowing to tak e t aking into a ccount of the multi-valuated arrivals are used-f&(= 
e xampl e. 

According to an implementation example, a zero-offset or a near-offset iso- 
15 offset collection is preferably constructed in stage -step a ). 

According to another implementation example, a fixed-velocity kinematic 
inversion technique such as a map migration is advantageously applied in stage -step 
c). 

According to another implementation example, the kinematic information is 
20 calculated in stage -step d) by tracing multi-offset rays on the interface of each 
interface-velocity pair. 
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According to another implementation example, sta§e-stee.f) is 
advantageously carried out by applying a prestack kinematic inversion method such 
as a prestack traveljime tomography. 

In its most general definition corresponding to the case where the number n of 
5 iterations to be carried out may possibly be zero, the method applies for example if 
the velocity distribution in each layer is hardly variable laterally and/or in the case 
where a sufficiently precise a priori on the velocity distribution in the layer considered 
is available. 

In i ts mor e r e strict i v e d e f i nition a more restrictive sense w here the number n 
10 of iterations to be carried out is at least 1 , the method applies for example in the 
case where the velocity distribution in the layers of the geologic formation is greatly 
variable laterally and/or in the case where no sufficient a priori on the velocity 
distribution in the layer considered is available. 

BRIEF DESCRIPTION OF THE F I GURES DRAWINGS 

15 Other features and advantages of the method according to the invention will 

be clear from reading the description hereafter of non limitative embodiment 
examples, with reference to the accompanying drawings wherein-: 

[[- ]]Figures 1a to 1c respectively show an application of the method to a common 
receiver type multi-offset collection (Fig..1a). to a common receiver collection on 
20 which are superposed the multi-offset traveljimes obtained with various layer 

velocity values (Fig._1b) and to a common receiver collection on which is superposed 
the predicted multi-offset traveljime curve showing the greatest coherence with the 
seismic event considered (Fig.J c). 
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DETAILED DESCRIPTION OF THE INVENTION 

Wg hav e se i sm i c Seismic records are obtained by means of a seismic device 
comprising a seismic source emitting seismic waves propagated in the subsoil, a set 
of seismic receivers are coupled with the medium which pick up the waves reflected 
5 by the subsoil discontinuities in response to the waves emitted, and of-a laboratory is 
used f or recording the picked up seismic signals p i ck e d up . 

I - Standard case 

W e pr e sent hereaft e r a A ,first implementation mode of the method Js 
described , in the case where the v e locit ie s v elocitv d istribution in the geologic 
10 formation preceding the layer i n qu es tion being examined is known or has been 

estimated othenA/ise, the geometry of the layer in question generating no triplications, 
the velocity to be determined being hafdiv -slightiv variable laterally or a sufficiently 
precise a priori on this velocity being available. The following stag e s steps a re 
carried out-: 

IS M From the prestack or multi-offset seismic records, con s truct i ng a collection 

of constant offset, is constructed preferablv zero-offset (in reality an approximation of 
a zero-offset collection, known as a^stack, characterized by a better signai-to-noise 
ratio than the individual iso-offset sections) or near-offset (corresponding to the 
lowest offset of the multi-offset records). 

20 If a coherent and sufficiently complete interpretation of the seismic event in 

progress can be carried out in this collection, picking Is directly performed in this 
collection. In the opposite case, in order to gain access to the traveUimes 
associated with the event in question, a detour can be made via the time or depth 
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migrated domain, followed by a time or depth demigration, as in the method 
described in-: 

[[- ]]Ehinger, A., and Lailly, P., 1995, Velocity mod e l det e rm i nat i on Model 
Determination by the SMART Method m e thod . Part 1-:. Theory-:. 65^^ Ann. Internat. 
5 Mtg., Soc. Expl. Geophys., Expanded Abstracts, pp.J39-742. 

1-2 A reference yelocity is selected for the layer in question, as well as an 
uncertainty oo-of this reference yelocity and a sampling interyal for the yelocity range 
thus formed. W e thu s Thus obta i n a series of yelocity samples for the layer in 
question is obtained : the reference yelocity can be any yelocity (that is. I iei-laterally 
10 and yertically yariable); in particular, if a yelocity distribution obtained othenvise is 
ayailable, this yelocity can be used as the reference yelocity; if no a priori information 
is ayailable for the yelocity of the layer, a laterally and vertically homogeneous 
reference velocity can , for example, be taken for exampl e. 

1-3 For each velocity sample, a fixed-velocity type kinematic inversion 
15 technique such as a map migration for example is applied in order to determine, from 
the traveUimes extracted from the constant-offset section, the geometry of said-the 
interface for the velocity sample in question. : w e thu s obta i n T hus a series of 
interface/velocity pairs for the event in question is obtained . A description of map 
migration can for example be found in the following reference-: 

20 [[- ]]Yilmaz, O., 2001 , Seismic Data Analysis - Processing. Inversion and 
Interpretation of Seismic Data-:. Society of Exploration Geophysicists. 

It may be impossible to find an interface model explaining the traveljimes for 
the fixed velocity in progress (impossibility of reaching convergence). .For a higher 
efficiency of the method, it is possible to select only interface/velocity pairs for which 



the difference between the travel Jimes observed and the traveljimes calculated on 
the model obtained after convergence is below a certain threshold. 

It can be noted that this fixed-velocity inversion can also be carried out by 
means of any other equivalent kinematic technique. 

5 1-4 The kinematic information associated with each interface/velocity pair 

obtained for example by means of multi-offset ray tracing on the interface of each 
interface/velocity pair thus obtained and for source-receiver pairs corresponding to 
multi-offset collections existing in the seismic records is calculated. ^Common 
midpoint collections are preferably used, in particular if the real velocity shows lateral 
10 variations that are not taken into account in the reference model. 

1-5 For each interface/velocity pair and for each multi-offset collection 
selected, the coherence between the multi-offset traveljimes thus calculated and the 
seismic records is evaluated.. A technique referred to as semblance calculation 
technique can for example be used, wherein the energy encountered in the seismic 

15 data along a curve that corresponds here to the multi-offset traveljimes calculated 
by ray tracing is summed. _ln an ideal case, there will be a velocity sample for which 
the predicted traveljime curve is perfectly superposed on the seismic event 
considered, which will translate into a maximum coherence measurement. Jt can be 
noted that this evaluation can be performed on a limited offset range, defined by the 

20 user or automatically, for example by examining the curve showing the evolution of 
the coherence as a function of the offset range. 

For each multi-offset collection, the travel Jime curve which best matches the 
seismic records is selected, and the multi-offset traveljimes forming this curve are 
reserved. .It can be noted that the travel.time curve exhibiting the best coherence 



can be adjusted if necessary by seeking in a predetermined vertical window the 
closest amplitude maximum (or minimum, depending on the polarity of the event 
considered) for each seismic trace of the multi-offset collection considered. 

1-6 The multi-offset travel.times thus collected for all the multi-offset 
5 collections selected are then injected into a prestack kinematic inversion method 
such as a traveljime tomography in order to determine the geometry and the velocity 
of the layer in question. 

I- 7 One may havo to i torato the The p revious stages mav have to be iterated if 
the coherence with the seismic records is not considered to be globally satisfactory 

10 by the operator. 

II - Particular cases 

II- 1 In cases where the velocity distribution in the layers of the geologic 
formation is greatly variable laterally and/or where no sufficiently precise a priori on 
the velocity distribution in the layer considered is available, the method according to 

15 the invention comprises the following steps stag e s : 

The previous stages are applied iteratively. the travel Jimes collected during 
an iteration being used for updating the velocity model by means of a prestack 
kinematic inversion method such as traveljime tomography, a model which is then 
used as the input velocity model for a new iteration of the method. 

20 It can be noted that, in order to be free from the effects of the lateral velocity 

variations that might have a smaller wavelength than the maximum offset of the 
multi-offset collections, the offset range considered by the method can be widened 
as the iterations progress.. More precisely, the method is initiated on a limited offset 



range, then the velocity distribution found from this limited offset range is used as the 
reference model for a new iteration of the method during which a wider offset range 
will be considered, and so forth,^ It can be noted that the offset range considered at a 
given iteration can be selected for example by examining the curve showing the 
5 evolution of the coherence as a function of the offset range. 

For the cases where the velocity distribution varies greatly laterally, it is 
possible to adopt furthermore a multi-grid approach, more precisely, the method is 
initiated on a loose grid of multi-offset collections (in order to determine the great 
wavelengths of the velocity variations), then multi-offset collection grids that are 
10 increasingly finer as iterations progress are taken into account. 

11-2 In the case where the sought interface geometry generates triplications 
and if it is also desired to determine the zones of the layer considered generating 
these triplications, it is first and foremost necessary to have inversion methods such 
as (fixed-velocity and prestack) kinematic inversion methods for example, allowing to 
IS take account of multi-valuated arrivals. Such a method is described in the following 
reference-: 

[[- ]]Delprat-Jannaud. F. and Lailly, P., (1995), How to Handleh afldle Multiple mu l t i p le 
Arrivals a rr i va ls Journal of Geophysical Research. 100. N0..BI, 703-715. 

The input data of this type of method are the multi-valuated travel.times 
20 associated with the seismic events considered, and the ray parameters associated 
with these traveUimes. In order to obtain the multi-valuated traveljimes and the ray 
parameters associated with the constant-offset data, a migration-demigration loop 
can be used (as described in the aforementioned method by Ehinger and Lailly 
(1995)) by carrying out the demigration stage by means of ray tracing allowing to 



c a lcu la t e calculation of multi-valuated arrivals. _Such ray tracing is for example 
described in-: 

[[- ]]Jurado, F., Lailly, P.. and Ehinger, A., (1998), Fast 3D two po i nt ravtrac i ng Two 
Point Ravtracinq fo r trav e lt i m e tomography T ravel Time Tomography : Proceedings 
S of SPIE, Mathematical Methods in Geophysical Imaging V, 3453, 70-81 . 

Then, after forming the various interface/velocity pairs by kinematic inversion 
at fixed velocity on the obtained t ime/rav parameters data thus obta i n e d , the 
kinematic information is calculated as described in the most general case by means, 
for example, of multi-offset ray tracing in each one of these models and for each 
10 multi-offset collection selected, but this time adapted to the multi-valuated arrival 
calculations as described in the aforementioned document by Jurado et al. (1998). 
Both prestack multi-valuated traveljimes and their associated ray parameters are 
thus obtained for each multi-offset collection selected. 

Then, the coherence analysis between the travel.time curves thus predicted 
15 and the events considered in the seismic data is carried out as described in the most 
general case, but independently for each branch of the multi-valuated arrivals.. The 
traveljimes thus collected for each branch of the multi-valuated arrivals are then 
injected, with their associated ray parameters, in a prestack kinematic inversion 
method allowing to tak e taking a ccount of the multi-valuated arrivals, as described in 
20 the aforementioned document by Delprat-Jannaud and Lailly (1995). 

W e hav e d e scr i b e d h e r e an A n application of the method to determination of a 
velocity model of seismic waves in an underground formation has been described . It 
is however obvious that the method can also apply to determination of the velocity of 



propagation of other types of waves in a heterogeneous model exhibiting 
discontinuities. 
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